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^ ' There has been a proof by Sewell that the hypothesis of off-diagonal long-range 

H : 

y3 , order in the reduced density matrix p2 implies the Meissner effect. We present in this 

note an elementary and straightforward proof that not only the Meissner effect but also 
the property of magnetic flux quantization follows from the hypothesis. It is explicitly 
shown that the two phenomena are closely related, and phase coherence is the origin 
for both. 
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Introduction 

Penrose and Onsager^ first proposed tliat tlie mathematical statement of Bose-Einstein 
condensation in tlie case of interacting boson systems is the existence of off-diagonal long- 
range order (ODLRO) in the reduced density matrix pi. Yang^ later made a thorough study 
of the subject and generalized the concept to interacting fermion systems. It is now generally 
accepted that the superconducting state of electrons is characterized by the existence of off- 
diagonal long-range order in the reduced density matrix p2- 

In view of the fundamental importance of the concept of ODLRO, it is important that 
there be a rigorous proof of the assertion that ODLRO in the reduced density matrix p2 
implies both the Meissner effect and the quantization of magnetic flux, which are the basic 
characteristic properties of superconducting states. There has been a recent proof by SewelF 
of the Meissner effect as a consequence of ODLRO, under certain simplifying assumptions. 
We shall extend in this note the proof to the case of magnetic flux quantization. 

The proof we present is simple and elementary, and clearly shows that both the Meissner 
effect and the flux quantization are closely related phenomena, and the underlying origin for 
both is the coherence of phases. 

System in Uniform Magnetic Field 

We consider a system of A^ electrons inside a crystal, which is placed in a uniform external 
magnetic fleld. We make the simplifying assumption that the magnetic fleld B inside the 
material is also uniform, with the corresponding vector potential, given by A(r), which is of 
the general form: 

A(r) = Ao(r) + V0(r), (1) 

where Ao(r) = |B x r, and the gauge 0(r) is, in general, multi-valued in r. The Hamiltonian 
for the system is 

H = E ^[^V, + ^A(r,)]2 + 1 E Vir,,) + ^ V;(r,), (2) 

J j,k j 

where Vc{rj) represents the interaction of the crystal lattice on the j-th electron, and the two- 
particle potential V{rjk), which includes the effective interaction due to phonon exchange. 



depends on the distance between the two electrons. Let 

^n(ri,r2,---,r7v), n = l,2,..., 

represent a complete orthonormal set of single-valued and properly antisymmetrized eigen- 
functions of the Hamiltonian 

H^pniri, r2, ■ ■ ■ , Tat) = ^^^^(ri, r2, ■ ■ ■ , tn). (3) 

We neglect the edge effects near the boundaries, and assume the usual periodic boundary 
conditions for the wave functions. The phenomenon of flux penetration near boundaries is 
beyond the scope of the present idealized considerations. 

Two Complete Sets of Eigenfunctions 

Now, consider a space displacement 

^j ~^ ^j ~ ^' 
under which 

A(r,) ^ A(r, - a) = A(r,) + V,[a ■ Ao(r,) + 0(r, - a) - 0(r,)]. (4) 

We observe that a space displacement induces a gauge transformation in the vector potential. 
It is this property that makes it possible for a simple proof of the flux quantization and the 
Meissner effect as consequences of the ODLRO hypothesis. The lattice potential Vc{rj) has 
the symmetry of the lattice, and is invariant under a displacement only if the displacement 
respects the symmetry of the crystal lattice. 

We flrst consider the case where the lattice potential Vc{rj) in the Hamiltonian is omitted, 
and will come back later to consider the case where it is not. Without the lattice potential, 
the eigen-equation (3), under the displacement Tj -^ Tj — a, becomes 

(E ^[-^^- + -M'^^ - ^)]' + 2 E yirjk)}Mri - a, r2 - a, ■ ■ ■ , rjv - a) 

= £^„^„(ri -a,r2-a, ■■■,rAf-a). (5) 



Noting that a displacement induces a gauge transformation, as in (4), in the vector potential, 
we can absorb the gauge into a phase for the electron wave function to obtain 

{Y. ^["^i + -^(^i)]^ ^ 2 ^ ^(^ifc)} ^^P[^ Y. Xa(rj-)]V^n(ri - a, rs - a, • ■ ■ , r^ - a) 

J j,k j 

= ^„exp[— ^Xa(rj)]^„(ri - a, ra - a, ■ ■ ■ , rjv - a), (6) 

where 

Xa(rj) = a ■ Ao(rj-) + 0(rj- - a) - 0(rj-)- (7) 

It is seen that '^/'^(ri, r2, ■ ■ ■ , r^r), given by 

<(ri, ra, ■ ■ ■ , r^v) = exp[— ^ Xa(rj)]^n(ri - a, rs - a, ■ ■ ■ , r^y - a), (8) 

are again eigenf unctions of the original Hamiltonian. In the following we will only con- 
sider the case of a to be infinitesimal on macroscopic scale^. In this case the phase Xa(r) 
is single- valued in r. Therefore, ipl^ form another complete orthonormal set of single- valued 
eigenfunctions. We thus have two bases of eigenfunctions, ipn and ipl^. 

Reduced Density Matrix for Two Particles 

The reduced density matrix element p2(i''i, i''2; ^i, r2) is defined in terms of a complete 
orthonormal set of energy eigenfunctions by 

P2(r'i, r'2; ri, rs) = y ■ ■ ■ y ["^ _ 2) f z ^ ^^M-En/kT) 

xipnir'i, r'2, rg, ■ ■ ■ , r7v)C(ri, ^2, rs, ■ " " , ^n), (9) 

where Z is the partition function, and the normalization is such that 

Trp2 =11 P2(ri, ^2] ri, r2)(iri(ir2 = N{N - 1). 

We note that P2(r'i, r'2; ^1, T'2) is single- valued, as the wave- functions are single- valued. Since 
the reduced density matrices are basis independent, we can equally represent 
P2(r'i, r'2; ri, r2) in terms of the primed set of energy eigenfunctions: 

P2(r'i,r'2;ri,r2) =]■■■] ("^T -2)\ Z^ exp(-K/A;T) 



xV^n(r'i, r'2, rs, ■ ■ ■ , r7v)^^*(ri, rs, rg, ■ ■ ■ , r^v), (10) 

Expressing the primed eigenfunctions in the above equation in terms of the original eigen- 
functions, we obtain 

P2(r'i, r'2; ri, r2) = exp{ — [xa(r'i) + Xa(r'2) - Xa(ri) - Xa(r2)]} 



xz/'„(r'i - a, r'2 - a, ra - a, ■ ■ ■ , Fat - a)?/'*(ri - a, r2 - a, ra - a, ■ • ■ , rAr - a). 



On account of the periodic boundary conditions satisfied by the wave functions, we can make 
an innocent shift of the integration variables in the integral in the above equation, resulting 
in the relation of the density matrix at different space points: 

P2(r'i,r'2;ri,r2) = exp{ — [xa(r'i) + Xa(r'2)-Xa(ri)-Xa(r2)]}p2(r'i-a, r'2-a; ri-a, r2-a). 

(11) 
This equation, as we shall see, is of basic importance. 

Off-Diagonal Long-Range Order 

We can represent the reduced density matrix element in the spectral form: 

P2(r'i,r'2;ri,r2) = E««'^«(r'i,r^)'5:(ri,r2), (12) 

s 

where as are eigenvalues, and $s(ri,r2) the properly normalized eigenfunctions of P2: 

<l'*(ri,r2)<l',(ri,r2)rfrirfr2 = 1. 
We are interested in the property of p2 in the off-diagonal long-range (ODLR) limit: 

|r'j-rfc|^oo, j,k = 1,2 

while keeping |ri — r2| and |r'i — r'2| finite. The hypothesis of ODLRO in p2 is the statement 
that in the ODLR limit the reduced density matrix element factorizes in the form 

P2(r'i, r'2; ri, r2) -^ ao^{r\, r'2)$*(ri, r2), (13) 



where ao is the largest eigenvalue of 0{N), and $ the corresponding eigenfunction^. 

Applying the factorization property of the ODLRO hypothesis to both sides of the relation 
(11), we obtain, in the ODLR limit, 

$(r'i,r'2)r(ri,r2) = 

exp{ — [xa(r'i) + Xa(r'2) - Xa(ri) - Xa(r2)]}$(r'i - a, r'2 - a)$*(ri - a, r2 - a), (14) 



which implies^ 



$(ri, r2) = /(a) exp{ — [xa(ri) + Xa(r2)]}$(ri - a, r2 - a), (15) 



where /(a) is a possible displacement dependent, but position independent, phase factor^, 
with /(O) = 1. Thus, the function $, which is to be identified with the Ginsburg-Landau 
wave function, at different space points is related by a phase factor. Repeating successively 
for infinitesimal displacements along a closed path C, and picking up products of phase 
factors and the path-dependent f-factors, we obtain from (15) the relation 

$(ri, r2) = F{C) exp{2^ / A(r) ■ rfr}$(ri, r2), (16) 

ch Jc 

where F{C) represents the product of the f-factors, and use has been made of the single- 
valuedness of the function $(ri,r2), a consequence of the single-valuedness of the density 
matrix. 

Meissner Effect 

We now use (15) to prove the Meissner effect. In a simply connected region we apply 
(15) to two sequences of successive displacements^, first a followed by b, and then b followed 
by a. Comparing the two resulting relations, we obtain 

is 
exp{ — [xa(ri) + Xb(ri - a) + Xa(r2) + Xb(r2 - a)]} = 
ch 



exp{^[Xb(ri) + Xa(ri - b) + Xb(r2) + Xa(r2 - a)]}. (17) 



Namely, 



2e / 
expU—rf A(r) • dr} = 1, 

mJ(axb) 



where the loop path of integration is a parallelogram formed by a and b. Therefore, 

ch 
B ■ (a X b) = n — , n = integers. (19) 

While a and b can be continuously changed, the right-hand side of (19), being proportional 
to an integer, can not. Consistency then requires 

B = 0. (20) 

This is the property of the Meissner effect. 

Flux Quantization 

Next, we use (15) and (16) to prove the quantization of flux, when the loop C is not 
simply-connected, like a path winding around a tunnel region inaccessible to the electrons. 

As we have already proved the Meissner effect, namely B = 0, it is clear from (16) that 
F{C) = 1 if C is simply-connected. To prove the quantization of flux, we need to prove 
F{C) = 1 for multiply-connected loops C. 

Now, the factor /(a) in (15) is independent of r. This implies that the property of F{C) 
is independent of the location of the loop C. Corresponding to a multiply-connected loop C, 
we can consider a loop of the same shape in a simply-connected region, for which F{C) = 1. 
It follows that F{C) = 1 also for C in a multiply-connected region. In general, therefore, 

F{C) = 1 (21) 

for any loop path C. (16) then becomes 

2e f 
$(ri, ra) = exp{i— j A(r) ■ rfr}$(ri, ra). (22) 

Consequently, there is the quantization condition 

f ch 

j> A(r) ■ dr = n — , n = integers, (23) 



The loop path C may or may not be simply-connected. When C in (23) is simply- 
connected, it can be altered infinitesimally and continuously, implying, again, the Meissner 



effect. When C is not simply-connected, lik;e a path winding around an inaccessible tunnel, 
the path can not be deformed across this inaccessible region. Then the quantization condition 
(23) implies that the total flux® passing through the loop C can only have quantized values. 
We note that in our idealized considerations the loop path C is away from the boundaries. 

We can prove F(C) = 1 more explicitly by observing some simple properties of /(a): 

1. As /(a) is independent of r, we can choose to explore the properties of /(a) by 
considering displacements in regions where B = 0. 

2. In any region of B = 0, consider the succession of two infinitesimal displacements a 
and b, as well as the displacement a + b. Because B = in the region enclosed by the three 
vectors a, b and a + b, it is clear from (15) that 

/(a + b) = /(a)/(b). (24) 

(24) offers another more explicit proof of (21). Recall that F{C) is the product of the 
f- factors along the loop path. Hence, for a closed loop, on account of (24), we again obtain 

F(C) = /(a)/(b)/(c)--- 

= /(a + b + c + ---) 
= /(£rfr) = /(0) = l. 



Inclusion of Lattice Potential 

The lattice potential V^(r), which is not invariant under arbitrary displacements, has so 
far been omitted in our consideration. We now consider the case when it is included. 

Let the symmetry of the lattice be defined by three basic lattice vectors: ei, e2, e3. The 
lattice potential Vc{r) is then invariant under the space displacement 

r, ^ r, - Ar, (25) 



if9 

Ar = riiBi + 71,262 + ^363, ni, n2, ris = integers. (26) 

Under such displacements, the total Hamiltonian is invariant, except for an induced gauge 
term in the vector potential, like that in (4). This is exactly like the situation we had 
when the lattice potential was omitted, except, of course, that the displacements have to be 
restricted to those in (26). All previous arguments apply, leading to the conclusion, similar 
to (19), that 

nh 

(27) 



(2^ 



B- (e^ X Gfc) =njk , 


rijk = integers 


which, in turn, implies that 




\/2e^'''^'^^' 


nj = integers 


where 




V = ei- (62 


X 63). 



Namely, the magnetic field B, if non- vanishing, inside the material can only be along a 
discrete set of directions. Now, following SewelF, we make the assumption that B is along 
the direction of the external applied field. As the direction of the external field can be 
arbitrarily chosen, the magnetic field B inside the material has to vanish, if not to contradict 
(28), giving rise to the Meissner effect. 

In the above proof we have made the assumption that the direction of B inside the 
material lies along the direction of the external field. We can, however, present another 
argument without having to make this assumption. Again, consider the requirement imposed 
by (28). Typical lattice dimension of superconducting materials is on the order of lA. The 
smallest value of B, if not zero, has to be on the order of 

B ~ {lA)-^ lO^G, (29) 

according to (28). This is many orders stronger than the typical critical field strengths 
of less than lO^G, beyond which the superconducting state no longer exists, for supercon- 
ductors. Therefore, to be consistent with (28), the only possible value of B is zero for a 
superconducting state. 

9 



The same reasoning that lead to the quantization condition (23) apphes, and leads to 

^A(r)-Ar = n— , n = integers, (30) 

c ^^ 

where loop path C consists of a succession of displacements of the type specified by (26). 
This quantization condition, when applied to multiply-connected paths, implies the phe- 
nomenon of flux quantization. 



Discussion 

We conclude with a few remarks. 

1. The function $(r,r) corresponds to the Ginsburg-Landau wave function^". The suc- 
cessive displacements taken in arriving at equation (16) clearly demonstrates the origin of 
the characteristics of the superconducting state: Coherence of phases is responsible for the 
properties of flux quantization and the Meissner effect. It is also made clear by the proof 
given in this note that these twin properties of superconductivity have a common origin, an 
understanding perhaps not sufficiently emphasized in usual presentations. 

2. The aspect in the hypothesis of ODLRO that has been made use of in the proof of fiux 
quantization and the Meissner effect is the factorization property^ of p2 as in (13), which 
is responsible for the coherence of phases in $(ri,r2). The ODLRO hypothesis implicitly 
assumes that in the sum in (12) there is only one term, having a large eigenvalue on the 
order of 0{N), that survives in the ODLR limit. 
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^G. L. Sewell, J. Stat. Phys. 61, 415(1990). 

^ In proving the flux quantization, we shall need to consider the case where there are 
regions inaccessible to the electrons, like a tunnel. In these inaccessible regions, the 
magnetic field does not necessarily vanish. In order to guarantee the single-valuedness 
of phase Xai^j) ^^ (7), we need to restrict to infinitesimal a, so that 0(rj — a) — 0(rj) 
is single-valued, although 0(rj — a) and (/){rj) individually are not. 

^Here $(ri,r2) has been renormalized, i.e., a factor ^^oq has been absorbed into it. 

^We thank W. Weisberger for pointing out the possible existence of this displacement- 
dependent, but r-independent phase factor. 

^The consideration parallels Sewell's, in Ref. 3 , in his derivation of the Meissner effect. 

^The total flux includes the flux through the tunnel as well as the flux of the penetrated 
magnetic field into the boundary layer surrounding the tunnel. 

^For reasons explained in footnote"^, we only consider Ar which is infinitesimal on macro- 
scopic scale. Thus, the integers rii, n2 and n^ can not be large. 

^"^It is the factorization property that was assumed by L.P. Gor'kov [JETP 7, 505(1958)] 
to derive his effective equations for superconductivity. 
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